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ABSTRACT 

Several lines of evidence suggest that the galaxy cluster C10024+17, an apparently relaxed system, 
is actually a collision of two clusters, the interaction occurring along our line of sight. Recent lensing 
observations suggest the presence of a ring-like dark matter structure, which has been interpreted as 
the result of such a collision. In this paper we present A^-body simulations of cluster collisions along 
the line of sight to investigate the detectability of such features. We use realistic dark matter density 
profiles as determined from cosmological simulations. Our simulations show a "shoulder" in the dark 
matter distribution after the collision, but no ring feature even when the initial particle velocity 
distribution is highly tangentially anisotropic {cre/cTr » 1). Only when the initial particle velocity 
distribution is circular do our simulations show such a feature. Even modestly anisotropic velocity 
distributions are inconsistent with the halo velocity distributions seen in cosmological simulations, and 
would require highly fine-tuned initial conditions. Our investigation leaves us without an explanation 
for the dark matter ring-like feature in CI 0024+17 suggested by lensing observations. 
Subject headings: cosmology: dark matter — galaxies: clusters: general — galaxies: clusters: individ- 
ual (CI 0024+17) — methods: N-body simulations 



1. INTRODUCTION 

Clusters of galaxies have proven to be interesting labo- 
ratories for studying the dynamics of the different kinds 
of matter in the universe. Early measurements of galaxy 
velocity dispersions in clusters demonstrated that most 
of the mass of galaxy clusters was in the form of a non- 
luminous component (|Zwickylll937l ). Modern cosmologi- 
cal constraints (e.g. Turner 2000) suggest that this non- 
luminous material must be mostly non-baryonic. In addi- 
tion. X-ray observations of clusters of galaxies reveal that 
most of the baryonic mass in clusters is in the form of a 
hot, diffuse. X-ray emitting gas, the intracluster medium 
(ICM). The presence of these different types of matter 
(stellar, gaseous, and dark) in these systems provides not 
only insights into the nature of the clusters themselves 
but also the nature of the different components of matter 
and their dynamics. 

Galaxy clusters are themselves the product of many 
mergers between galaxies, galaxy groups, and smaller 
clusters of galaxies (Davis et al. 1985). Observations in- 
dicate that this process is still ongoing in many systems. 
A famous example is the so-called "Bullet Cluster" (IE 
0657-56), a merging system in which X-ray and weak 
lensing observations demonstrate a clear separation be- 
tween a collisionless dark matter component and the in- 
tracluster medium (|Markevitch et al.l 120021 : 1 Clowe et al.l 
[2006). 

A system that has also garnered attention recently 
is CI 0024+17, a cluster at z = 0.395 with weak and 
strong lensing. Early attempts at reconstructing the 
mass profile of thi s system using lensing (Tvson et al. 
119981 : iBroadhurst e t al. 2000; Come rford et al. 2006) sug- 
gested a conflict with the predictions of the standard 
cold dark matter (CDM) model due to the flatten- 



ing of the density proflle in the inner regions of the 
cluster. Cosmological simulations assuming CDM in- 
dicate that galaxy cluster mass proflles should exhibit 
a logarithmic slope in the inner regions (e.g. Navarro, 
Frenk, & White 1997). This apparent discrepancy 
(and others) led to suggestions that the CDM paradigni 
would need to be modified (SDergel & Steinhardt 20Q0|; 
Hogan Dalca nton 20001: IMoore et al.ll200a) . for exam- 
p le to include self-inte raction of the dark matter. 

ICzoske et all (|200lL |2002) demonstrated that the red- 
shift distribution of the cluster galaxies in CI 0024+17 
is bimodal and suggested that the system is composed 
of two clusters undergoing a collision along the line of 
sight. They also performed a simulation demonstrating 
that such a scenario reproduces not only the bimodal 
redshift distribution but also the fiattening of the central 
density profile. Observations of the cluster by Chandra 
(Ota et al. 2004) and XMM-Newton (Zhang et al. 2005) 
revealed that the surface brightness profile is better fit 
by a superposition of two ICM models rather than one. 
They suggested on the basis of the isothermal tempera- 
ture profile that the system had returned to equilibrium 
after the collision and that consequently the encounter 
must have occured several Gyr ago. 

Recentl y, a^ weak lensing analysis presented by 

iJee et all (|2007 ) revealed a ring-like structure in the pro- 
jected matter distribution. They proposed that dark 
matter from the cores of the clusters had been disrupted 
and ejected from the systems by the collision, and they 
demonstrated that such features could be reproduced in 
a simulation of a collision of two pure dark matter halos. 
On the basis of this result they suggested the current 
state of the system is ~l-2 Gyr past the pericentric pas- 
sage of the cluster cores. 

We wish to shed additional light on the formation 
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of ring structures within the dark matter distribution. 
"Ring galaxies" provide another astrophysical situation 
in which "particle" rings are produced due to collisions 
(e.g., the Cartwheel Galaxy). We look to the previous 
A^-body simulations of these phenomena to guide our un- 
derstanding of what might lead to such a feature in the 
dark matter in clusters of galaxies. Two major differ- 
ences exist between the dark matter particles in galaxy 
clusters and star "particles" in disk galaxies. In the case 
of disk galaxies, the stars are concentrated in a disk, and 
the velocities are essentially circular. In galaxy clusters, 
dark matter particles form a roughly spheroidal distri- 
bution and have isotropic to radially anisotropic velocity 
dispersions. To investigate the effects of the second of 
these differences we perform a set of pure dark matter 
(A^-body only) simulations with a physically motivated 
dark matter profile and variations on the velocity distri- 
bution of the dark matter particles. 

Throughout this paper we assume a flat ACDM cos- 
mology with h — 0.7 and = 0.3. 

2. SIMULATIONS 
2.1. Method 

We performed our simulations using FLASH, a parallel 
hydrodynamics/ A^-body astrophysical simulation code 
developed at the Center for Astrophysic al Thermonu- 
clear Flashes at the University of Chicago (jFrvxell et al.l 
[2000). FLASH uses adaptive mesh refinement (AMR), a 
technique that places higher resolution elements of the 
grid only where they are needed. AMR allows us to 
achieve high resolution these regions without having to 
fully resolve the whole grid. FLASH includes an TV-body 
module which uses the particle-mesh method to solve 
for the forces on gravitating particles. The gravitational 
potential is c omputed usin g; a multigrid solver included 
with FLASH (|Rickedl2008l ). 

2.2. Initial Conditions 

For ease in setting up the par ticle distribution func- 
tions, we use a Hernquist profile (|Hernquistl 1 1 990l ) : 



Pdm(^) Ps 



(1) 



' r/a(l + r/a)^ ' 

where ps = Mo/(27ra^) is the scale density of the pro- 
file. The Hernquist mass profile that converges to Mq as 
r ^ oo: 



M(r) = Mo- 



'(r + a)2 

This form of the dark matter density profile is cho- 
sen because of the mathematical simplicity of the corre- 
sponding distribution functions (as will be shown below) 
and its resemblance to the NFW profile in that as r ^ 0, 
p{r) (X r~^. 

To determine the effects of a varying velocity 
anisotropy on the dark matter features in the simula- 
tion, we allow the 3D velocity dispersion of the parti- 
cles to vary from the isotropic form ar = cfq = ac/y. We 
assume ag = (Jcf) and parametrize the anisotropy using 

= 1 — which is taken to be constant over the en- 
tire radial range. Table [1] shows the values of (5 that were 
used for each simulation, and the corresponding ratio — . 



TABLE 1 
Velocity Anisotropy 



Simulation 




(To/ar 


E (erg) 


SI 


1/2 


0.707 


-4.45 X 10^3 


S2 





1 


-4.38 X 10^3 


S3 


-3/2 


1.58 


-4.08 X 10^3 


S4 


-5/2 


1.87 


-4.98 X 10^3 


S5 


-17/2 


3.08 


-5.28 X 10^3 


S6 


— oo 


oo 


-4.88 X 10^3 



In order to initialize the particle velocities, we sample 
the particle distribution function (hereafter DF) directly. 
The DF ^(r, v) is assumed to obey the following relation: 



pir) 



v)(iv 



(3) 



The DF of any steady-state, spherically symmetric 
system has a dependence on the phase space coordi- 
nates that comes in only through the "integrals of mo- 
tion" £ and L, where 8 = — \v'^ is the relative en- 
ergy and L = rvt is the angular momentum of a par- 
ticle (|Binnev fc Tremaind 11987). In these expressions, 
?/^(r) = —4){r) is the relative gravitational potential, and 
Vt is the tangential velocity. 



TABLE 2 
Initial Cluster Parameters 



Cluster 


Mo (Mo) 


R (kpc) 


a (kpc) 


Ps (MQkpc-^) 


1 


6.0 X 10^^ 


2000.0 


400.0 


1.49 X 10^ 


2 


3.0 X 10^^ 


1000.0 


200.0 


5.97 X 10^ 



For a constant velocity anisotropy, the DF takes the 
specific form (Binney & Tremaine 1987): 



(4) 



For a He rnquist profile, the ene rgy-dependent part of 
the DF is (iBaes fc Deionghdl2002D 



fs{£) = 



2/3 



r(5 - 2/3) 



(27r)5/2 r(l-/3)r(|-/3) 
X2F1 [5-213,1 -2(3,1 -I3;£ 



£.5/2-/3 



(5) 



where 2F1 is the hypergeometric function. For half- 
integer values of /3 this function can be expressed in terms 
of rational functions. To ease our investigation of this pa- 
rameter space we therefore choose values of P from this 
set, specifically the values (3 = 1/2, -3/2, -5/2, -17/2. We 
also include an isotropic setup (p = 0) and a setup where 
all velocities are initially circular {/3 = —00). Given the 
mass density function and the distribution function, we 
derive initial positions and velocities for the particles us- 
ing the acceptance-r ejection meth od. 

Finally, following ICzoske et al.l (|2QQ2l ) and 'Jee e t al] 
(|2QQ7l ) we assume a mass ratio of 2:1 for the clusters. 
In all of our simulations the clusters are initialized so 
that their centers are separated by the sum of their re- 
spective radii R (3 Mpc), and they are given an initial 
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relative speed '^rei = 3000 km/s, which is the inferred rel- 
ative velocity of the two co mponents of C10024 +17 seen 
in the redshift histograms of Czosk e et alJ (|2002i) . In that 
study they showed that a initial velocity of ^3000 km/s 
in their collision simulation reproduced this observed dis- 
tribution of redshifts. This value for the velocity is ap- 
proximately 2vff for the clusters, where Vff is the free- 
fall velocity from infinity. Such an initially high relative 
velocity n iight be difficult to achieve in a ACDM universe 
(see, e.g. iHavashi fc White! 12006). However, our choice 
is motivated by our desire to match the relative inferred 
velocity of the cluster components of CI 0024+17. The 
values of the halo parameters are given in Table O 

We refine the adaptive mesh using the dark matter 
density to resolve the cores of the clusters and other over- 
densities. For our box size of 10h~^ Mpc we achieve a 
minimum zone spacing of Ax = 9.77h~^ kpc. 

3. RESULTS 

For our set of simulations we have used the same cluster 
profiles, masses, initial separation, and relative velocity, 
but we have varied the anisotropy of the dark matter 
velocity dispersion. In each simulation, from t = 0.0- 
0.9 Gyr the dark matter cores accelerate toward each 
other and pass through each other. The resulting sudden 
impulse on the dark matter particles boosts their energies 
and causes a significant fraction of the mass in the cores 
to be ejected from the centers. Figure [1] shows slices 
through the z = coordinate plane of the dark matter 
density in the /3 = —5/2 simulation at different epochs. 
After the pericentric passage the outer regions of both 
halos have expanded, with a trail of dark matter strung 
between the centers. It can be seen that the smaller 
halo is more disrupted, with a larger amount of mass 
redistributed to the outer portion of the halo. 

Figures [2] through [3 show slices though the z = co- 
ordinate plane of the dark matter density, the projected 
dark matter density, and the projected density profile for 
each simulation at the epoch t = 1.5 Gyr 0.6 Gyr af- 
ter the collision, defined as the time when the centers of 
the halos are coincident). At this epoch, our simulations 
show no ring feature (defined as a "bump" in the radial 
distribution of dark matter particles projected onto the 
sky; i.e. an increase followed by a decrease in this distri- 
bution) for initial velocity distributions that are radially 
anisotropic {/3 = 0.5) and isotropic orbits {/3 = 0). For 
initial velocity distributions that are increasing tangen- 
tially anisotropic {/3 = —0.5, —2.5, —8.5), our simulations 
show a more pronounced shoulder but no ring. Only for 
an initially circular velocity distribution {f3 = — oo) do 
our simulations show a ring. Figure [8] shows the angle- 
averaged projected density profile at t = 1.5 Gyr for all 
of the different initial values of p. The initial profile is 
also shown for comparison. 

Even in the simulation in which our simulations show 
a ring (/? = — oo), the feature is short-lived. Figure [8] 
shows that the feature is prominent at epochs t = 1.5 
Gyr and t = 2.0 Gyr 0.6 and ^ 0.95 Gyr after the 
collision). However, hy t = 2.5 Gyr and t = 3.0 Gyr 
(^ 1.6 and ^2.1 Gyr after the collision), it has become 
insignificant. 

4. DISCUSSION 



4.1. Relevance to Observations of CI 0024+17 

Our simulations show that for physically and obser- 
vationally motivated dark matter particle distributions 
(NFW-like) and initial dark matter velocity distributions 
with a wide range of angular an isotropies , nodng fea- 
ture forms as in the scenario of iJee etldl (|2007[ ). The 
collision scenario for CI 0024- hl7, as mentioned b efore, 
is supported by the results of ICzoske et al.l (2002), who 
demonstrated the existence of a bimodal redshift distri- 
bution in the cluster galaxies. They also showed that a 
simulation of such a collision can explain the observed 
fiattening of the density profile in the cluster core even 
if the initial profiles are cuspy. However, they did not 
mention the existence of a ring feature in their paper. 

iJee et all (|2007[ ) presented results of a lensing anal- 
ysis of CI 0024-hl7. In addition, they presented the re- 
sults of simulations of head-on collisions of galaxy cluster 
dark matter halos. The dark matter density distributions 
of the galaxy clusters in their study were not NFW-like 
with central cusps, but softened isothermal spheres with 
p{r) (X [1 + (r/rc)^]~^; the initial velocity distribution of 
the particles was not reported. In these simulations, a 
clear ring-like feature is evident. 

In an attempt to understand the diffe rence between th e 
results of our simulations and those of iJee et al.l (|2007[ ), 
we have performed a simulation assuming a softened 
isothermal sphere for the density profile with the same 
parameters, and the same cluster masses, centr al densi- 
ties, a nd initial relative velocity as quoted in Jee et all 
(|2007f ). Unfortunately, jee et al 07 do not state the 
angular velocity distribution they use. We assume an 
isotropic velocity distribution for the dark matter parti- 
cles, in agreement with the velocity distribution in the 
inner regions of dark matter halos seen in ACDM sim- 
ulations of structure formation. Similarly to the /3 = 
case with the Hernquist profile, we find no evidence 
for a shoulder, let alone a ring-like feature (see[9|). This 
shows, as one might expect, that a modest difference in 
the density profiles between the Hernquist profile and the 
softened isothermal sphere profile is not the origin of the 
difference between the two results. 

Our simulations show no ring-like feature even for ini- 
tial particle velocity distributions that are highly tan- 
gent ially anisotropic; only for an initially circular ve- 
locity distribution do our simulations show a ring. As 
we discuss below, even modestly tangentially anisotropic 
velocity distributions are not expected for the dark mat- 
ter particles in galaxy clusters. This, together with the 
results of our simulations imply that the production of 
ring-like features in galaxy cluster collisions is likely to 
be rare. 

4.2. Comparison with ^'Ring Galaxies 

The phenomenon of "ring galaxies" (e.g., the 
Cartwheel Galaxy) is explained by a similar mechanism, 
that of a smal ler, more compa ct galaxy colliding with 
a disk galaxy (^Lvnds fc Toomrea J976) . A/'-body simula- 
tions of galaxy-galaxy collisions have demonstrated this 
phenomen on, for which the collisionless " particles" are 
stars (Thev s fc SDiegellll977tlToomrelll978[ ). In such sim- 
ulations, the velocity distributions of the stars are set up 
to be tangentially anisotropic, or < 0, and most of the 
orbits are circular. Also, the stars are concentrated into a 



Fig. 1. — Slices through the z = coordinate plane of dark matter density at the epochs t = 0.0, 0.5, 1.0, 1.5, 2.0, 2.5, and 3.0 Gyr for 
the f3 = —5/2 simulation. Each panel is 10 Mpc on a side. 





disk-shaped struc t ure rather than a spheroidal structure. 
iLvnds fc Toomrd (|197 6) demonstrated that such a colli- 
sion causes a crowding of particles on tangential orbits. 
In our simulations, as P is made more negative a more 
pronunced shoulder appears in the projected dark matter 
density, but a ring appears only when the initial velocity 
distribution is circular. However, ACDM simulations of 
structure formation have demonstrated that in the inner 
regions of dark matter halos the particle velocities are 



nearly isotropic {/3 ~ 0), and they become more radially 
aniso tropic in the outskirts of the halo (|Cole fc LacevI 
Il996l ). Our results show that under such conditions no 
ring forms. If the formation of dark- matter ring features 
in collisions between galaxy clusters is dependent upon 
tangentially anisotropic velocities of the particles, then 
it would require fine-tuned, high-angular momentum ini- 
tial conditions in the dark matter distribution, which is 
unrealistic under the standard CDM paradigm. 
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Fig. 4. — t = 1.5 Gyr snapshot of the [3 = -3/2 simulation. Left: Shce through the z = coordinate plane of dark matter density. Center: 
Dark matter density projected along the collision axis. Right: Projected density profile. 




-3.0 -2.0 





Fig. 5. — t = 1.5 Gyr snapshot of the /3 = -5/2 simulation. Left: Slice through the z = coordinate plane of dark matter density. Center: 
Dark matter density projected along the collision axis. Right: Projected density profile. 
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Fig. 6. — t = 1.5 Gyr snapshot of the f3 = -17/2 simulation. Left: Slice through the z = coordinate plane of dark matter density. 
Center: Dark matter density projected along the collision axis. Right: Projected densityprofile. 





1.0 2.0 3.0 




-3.0 -2.0 
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Fig. 7. — t = 1.5 Gyr snapshot of the (3 = — oo simulation. Left: Slice through the z = coordinate plane of dark matter density. Center: 
Dark matter density projected along the collision axis. Right: Projected density profile. 
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Fig. 8. — Left: Angle- aver aged DM density profiles at t = 1.5 Gyr into the simulation (~0.6 Gyr post-collision) for varying [3. As [3 
is made more negative a more pronounced shoulder-like feature appears in the density profile but a ring does not appear except in the 
(3 = — oo case. Right: Angle- aver aged DM density profiles at the epochs t = 1.5, 2.0, 2.5, and 3.0 Gyr for the (3 = — oo simulation. Note 
that the ringlike feature appears prominently at t = 1.5 Gyr (~0.6 Gyr after the collision) but is a transient feature that has completely 
disappeared at later times. 




Fig. 9. — t = 1.5 Gyr snapshot of a cluster collision with softened isothermal sphere density profiles and [3 = 0. Left: Slice through the 
z = coordinate plane of dark matter density. Center: Dark matter density projected along the collision axis. Right: Projected density 
profile. 
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5. CONCLUSIONS 

A recent lensing analysis of the cluster CI 0024+17 re- 
vealed possible evidence of a ring-like dark matter struc- 
ture. This has been interpreted as the result of a collision 
between two clusters of galaxies along the line of sight, 
a scenario for which there are other lines of evidence in 
this system. We have performed simulations of collisions 
between galaxy cluster dark matter halos to test this hy- 
pothesis, using density profiles for the dark matter halos 
motivated by observations and simulations and investi- 
gating a parameter space of initial velocity distribution 
where we vary the initial velocity anisotropy. Our sim- 
ulations show that a more pronounced "shoulder" fea- 
ture appears in the projected dark matter density when 
we make the velocities of the dark matter particles more 
tangentially anisotropic, by analogy with the phenomena 
of "ring galaxies." Our simulations show that, although 
the collision ejects a large amount of mass from the dark 
matter cores, a ring-like features does not form, even for 
initial velocity distributions that are highly tangentially 
anisotropic. Only when the initial velocity distribution 



is circular does a ring form. Since we do not expect 
dark matter particles in clusters of galaxies to have tan- 
gentially anisotropic velocities, our investigation of this 
parameter space leaves us without an explanation for the 
ring-like feature that appears in the dark matter distri- 
bution of CI 0024-hl7. 
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